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Abstract 

We have investigated the X-ray scattering signal of highly aligned multilayers 
of the zwitterionic lipid l,2-dimyristoyl-sn-glycero-3-phosphatidylcholine containing 
pores formed by the antimicrobial peptide alamethicin as a function of the pep¬ 
tide/lipid ratio. We are able to obtain information on the structure factor of the 
pore fluid which then yields the interaction potential between pores in the plane of 
the bilayers. Aside from a hard core with a radius corresponding to the geometric 
radius of the pore, we find a repulsive lipid-mediated interaction with a range of 
about 30 A and a contact value of 2.4 ksT. This result is in qualitative agreement 
with recent theoretical models. 
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1 Introduction 


The elucidation of lipid-mediated interaction forces between membrane proteins and the 
corresponding lateral distribution in the plane of the membrane is an important step 
towards a quantitative understanding of the functional mechanisms of membrane pro¬ 
teins and membrane peptides. Experimentally, the lateral structure and organization of 
multi-component membranes is as important as it is difficult to probe. While fluores¬ 
cence microscopy in biological and model membranes can be used to monitor domains 
and the distribution of proteins typically at the micron scale ([H) and down to a few hun¬ 
dred nanometers at best, atomic force microscopy can resolve lateral structures down to 
a nanometer (j^, but only in relatively stiff systems and rarely in the fluid state of the 
membrane. To this end, x-ray or neutron scattering from aligned fluid bilayers is an ex¬ 
cellent tool to probe correlations between proteins or peptides in the bilayer. In this work 
we show how the lateral and vertical intensity profiles of a peptide pore correlation peak 
can be analyzed as a function of peptide concentration to determine the corresponding 
interaction forces. 

A well-known example of biological function deriving from lipid-peptide interaction 
and self-assembly is the activity of a family of short and amphiphilic membrane active 
polypeptides denoted as antimicrobial peptides. These molecules bind to microbial cell 
membranes, subsequently causing an increase in membrane permeability and cell lysis. 
One such molecule is alamethicin, a 20 amino acid peptide from the fungus Trichoderma 
viride; it is well known that alamethicin acts by creating pores in the cell membrane (j^. 


Dv a convergence of multiple experimental techniques: 
(j3, lEI) and nuclear magnetic resonance (NMR) (j^ have 


This conclusion has been reached 
oriented circular dichroism (OCD) 
shown that, above a certain concentration, the binding state of the peptides changes from 
adsorbed parallel to the membrane to inserted into the bilayer. At the same time, a 
notable increase in membrane conductivity (0, and references therein) and permeability 


was measured. Pore formation is usually a highly cooperative process (jSl, 110|); this was 
conhrmed for alamethicin and a membrane-mediated interaction between peptides was 
invoked to explain the phenomenon (0). 

Although determining the interaction between (adsorbed or inserted) monomers is very 
difficult and, to our knowledge, has never been attempted, the interaction between already 
formed pores within the membrane can be studied using neutron or X-ray scatte ring from 
oriented multilamellar stacks, a method pioneered by Huang and collaborators f)lll- ll3l) . 

For the case of alamethicin, they observed a lateral correlation peak, which was at¬ 
tributed to liquid-like ordering of pores in the plane of the membrane and was modelled 
based on hard disk interaction, with very satisfactory results. However, in theses studies 
at most two peptide-to-lipid concentrations P/L were investigated for each system. 

Building upon this work, we gathered detailed information on the quasi two-dimensional 
fluid of pores in the lipid bilayer, using high-resolution synchrotron scattering from aligned 
multilamellar stacks of alamethicin/DMPC mixtures. We measured the two-dimensional 
scattering distribution for an entire concentration series P/L and performed a simultane¬ 
ous lineshape analysis on all recorded curves. 

We found that the in-plane interaction potential consists of a hard core, with a radius 
that agrees very well with the geometrical outer radius of the pore, and an additional 
repulsive contribution which can be described as a Gaussian, with a range of 31.5 A and a 
contact value of 2.41 ksT. The results are in qualitative agreement with recent theoretical 
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models (jl4 IlSf) . 

In principle, this method is readily applied to any peptide/lipid system, provided that 
well-oriented multilayer samples can be prepared. Thus the role of different parameters 
such as bilayer composition, temperature, nature of the aqueous medium etc. can be 
systematically studied. 


2 Materials and Methods 


2.1 Sample preparation and environment 

The lipid l,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC) was purchased from 
Avanti Polar Lipids Inc. (Birmingham, AL) with a purity of at least 99% and alamethicin 
was bought from Sigma Aldrich with a purity of at least 98.9%. Without further purifica¬ 
tion, the products were dissolved in a TFE/CHCI 3 , 1:1 vol/vol mixture at a concentration 
of 60 mg/ml for the lipid and 15 mg/ml for the peptide. The stock solutions were then 
mixed (and solvent added as necessary) to give the desired molar lipid/peptide ratio P/L, 
at a final lipid concentration of 20 mg/ml. The resulting solutions were then kept at 4°C 
for at least 24 hours before preparing the saniples. More details on sample preparation 
and on the choice of solvents are given by fjl 6 l-ll 8 l). 

Rectangular silicon substrates (15 x 25 mm^) were cut from 0.4 mm thick commercial 
wafers (Silchem Gmbh, Freiberg, Germany) and cleaned by sonicating them (during 15 
min) twice in methanol and then twice in ultrapure water (specific resistance > 18 MG cm, 
Millipore, Bedford, MA). Finally, they were abundantly rinsed in ultrapure water and 
dried under nitrogen flow. 

An amount of 0.2 ml of the solution was pipetted onto the substrates under a laminar 
flux hood, where they were subsequently left to dry at room temperature for a few hours 
and then exposed to high vacuum at 40 °G overnight to remove any remaining traces of 
solvent. They were finally stored at 4°G until the measurement. From the amount of 
lipid deposited, the thickness of the film can be estimated at about 3000 lipid bilayers. 

Before the measurement, the samples were placed in the experimental chamber and 
the hydrating solution was gently injected so as to avoid washing the lipid film off the 
substrate. Two types of sample chambers were used, the first one machined out of plexiglas 
and with an optical path of about 1.7 cm, and the second one made of teflon and with an 
optical path of about 1.1 cm. Both chambers have 0.3 mm thick kapton windows and were 
mounted on a metal heating stage temperature-controlled by water flow from a heating 
bath (Julabo Gmbh, Seelbach, Germany). 

For all P/L values, the hydrating solution was 100 mM NaGl brine containing 31 % 
w/w PEG 20000 (Fluka Ghemie Gmbh, Buch, Switzerland), yielding an osmotic pressure 
of approximately 1.68 10 ® PeQ. Additionally, for P/L = 1/12.5 we also performed mea¬ 
surements at 12.1 and 5.8 % PEG concentration, corresponding to 1.2 10® and 3.5 10^ Pa, 
respectively. The temperature was kept at 30 °G for all experiments. 


^This value was obtained from the data of Prof. Peter Rand, at the Membrane Biophysics Laboratory 
of Brock University, Canada: http://aqueous.labs.brocku.ca/osfile.html 
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2.2 Measurement 


The measurements were performed at the undulator beamline IDl of the European Syn¬ 
chrotron Radiation Facility (ESRF, Grenoble, France). The photon energy was set at 
19 keV by a double-bounce Si(lll) monochromator and the higher-order harmonics were 
cut by reflection on two Rh-coated mirrors. At this energy, the transmission of 1 cm of 
water is of 0.45, so the presence of the experimental cell does not pose any attenuation 
problems. 

Three types of measurements were performed: a) reflectivity scans (in the vertical 
scattering plane) up to a z component of the scattering vector ~ 0.8 A (see Fig . 


4]) give access to the electronic density profile of the bilayers along the z direction (jl9 
20l ). However, they correspond to averaging over the plane of the bilayer, so all lateral 
information is lost, b) CCD images are taken (using a Peltier-cooled camera, 1242 x 1152 
pixels, from Princeton Scientific Instruments Inc., New Jersey, USA) at a hxed incidence 
angle of the X-ray beam onto the sample and correspond to sections through the 
reciprocal space with the Ewald sphere (see Fig. E]); they provide a global image of the 
q-space and the position of the pore signal can be quickly determined, c) Quantitative 
measurements were performed using a point-detector (Cyberstar scintillation detector 
from Oxford Danfysik, Oxford, UK). Transversal (along qy) scans were taken through the 
pore scattering signal, with wide open slits in the vertical direction, covering a q^ range 
between 0.14 and 0.18 A . For some samples, longitudinal scans (along g^) were also 
taken. Their trajectories in g-space are shown in Figure [3] (right) as dotted lines. 


2.3 Analysis 


“pure” 


The alamethicin pores are dispersed in the lamellar phase matrix. Since the 
lamellar phase gives a signal confined to the vicinity of the Bragg peaks, from the Babinet 
principle it ensues that the off-axis scattering is the same as for a system where the 
density prohle of the lamellar phase is subtracted, and one is left with hctitious “pore - 
bilayer” objects in a completely transparent medium. Furthermore, as the pores represent 
a collection of identical and similarly oriented objects (up to an azimuthal averaging), the 
classical separation of the scattering intensity in a structure factor multiplied by a form 


factor can be applied (1211 1. yielding: J(q) = S'(q) ■ |F(q)| , with: 


S{qz,qr) = ^ 


N-l 


exp (-iqrfc) 


k=l 


( 1 ) 


where N is the number of objects and object “0” is taken as the origin of the coordinates. 
If there is no in-plane ordering, S only depends on the absolute value of the in-plane 
scattering vector g^ = ^g^ -|- g^. The form factor is given by: 


i-d/2 


r-R 


F{qz,qr) = Tj / dze / drrJo(grr) / d6p{r,6,z) 


V 


( 2 ) 


'-d/2 


where p is the electron density and V is the integration volume (correponding to the size 
of the object). 

The first step in computing the structure factor is determining the numerical density of 
pores in the plane of the bilayer (or, conversely, the area per pore). For the alamethicin, we 
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use the values given in the literature for an 8-monomer pore in DLPC fjl2l) : the peptides 
are modelled as cylinders of 11A in diameter; however, the effective cross-section of a 
peptide is only 66 (the rest being occupied by lipid chains). For the DMPC (j^ Table 
6 ), we consider that the area per molecule at 30 °C and in the absence of applied pressure 
is Ao = 59.6 A ; in our experiments, the osmotic pressure reduces it to A = 59 A , the 
area compressibility modulus being Ka = 0.234 N/m fli^. Table 1). We only detect a 
signihcant off-axis signal (assignable to the pores) for P/L > 1/25, concentration at 
which more than 80% of the peptide is in the inserted state, and this ratio increases to 
more than 90% for P/L = 1/20 (js], |2^ [^; thus we consider that all alamethicin molecules 
are involved in pore formation. 

The area per pore is determined assuming that the peptides are straight cylinders 
placed at the vertices of a regular hexa-, hepta- or octagon. The water pore is the 
inscribed cylinder (tangent to the peptides), with a radius Rim = 5.5, 7.2 and OA for 
m = 6, 7 and 8, respectively: 


A — 

^pore 'I I'irn, 


-|- m X 66A^ -|- 


m 1 


X 


59A^ 


(3) 


2 P/L 

where the 1/2 factor in the third term (corresponding to the area taken by the lipid 
molecules) accounts for the two monolayers. For simplicity, the kink in the peptide (see 
section ED and the polydispersity in aggregation number (12511 were neglected. 

To determine the form factor, we used the molecular dynamics (MD) result^ of Tiele- 
man et al. (1261) . who studied alamethicin pores of different sizes in a POPC bilayer. The 
form factor was computed according to Equation ([2]), for a 30 x 30 A patch containing 
the pore and for a similar patch containing only lipids (obtained by tiling four times a 
15 X 15 A^ patch from the same simulationjf]. The effective form factor used is the dif¬ 
ference of the two. We neglected the difference between POPC and DMPC when using 
the resulting form factor in our hts. Since the MD simulations indicate that the hexamer 


is the most stable conhguration in POPC (l26fl and neutron scattering results hnd 8-9 
monomers per pore in DLPC (jl2l) . we considered the pore conhgurations with 6, 7 and 8 
monomers. As we shall see later, the 7-monomer conhguration gives the best hts, so all 
results presented in the following correspond to this conhguration. 

In Figured] we present a (y,z) cross-section through the form factor |T(q)| of a pore, 
after subtraction of the pure bilayer background and azimuthal averaging. Directions x 
and y are equivalent. 

For a visual representation of the scattering object (pore - bilayer) we performed a 
Fourier transform of T(q) back to real space, yielding the density prohle shown in Figure 
dJ The peptide monomers are clearly visible as higher density streaks. 


^These results were obtained for DPhPC; The peptide is in the inserted state at all measured concen¬ 
trations in DLPC (12). 

^Available on the web site of Dr. Peter Tieleman, Dept, of Biological Sciences, University of Calgary: 
http://moose.bio.ucalgary.ca/Downloads/ (used with permission). 

"‘The patch sizes are 40 x 40 A/ and 20 x 20 A , respectively in the case of the 8-monomer pore. 
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Results 


3.1 Structure of the scattered signal 

To serve as an illustration for the discussion of the results, we show in Figure |3] a diagram 
of the reciprocal space structure for a multilamellar stack on solid support, as well as an 
actual CCD image (which amounts to a cut by the Ewald sphere). The off-axis signal 
(exhibiting a maximum around = 0.1 A ) is due to the alamethicin pores; for clarity, 
it is not represented in the diagram on the left. In order to bring up this (very weak) 
signal, the Bragg sheets are severely overexposed. The image was taken at an incidence 
angle a* = 0.55°, for a sample with P/L=l/20. 


3.2 Perfectly aligned samples 


We checked the quality of the samples and their alignment by performing reflectivity 
measurements fllOl . l20l ) . The mosaicity can be estimated at about 0.01° from the rocking 
scans. The reflectivity curves are shown in Figure H] for four different P/L values. Seven 
Bragg peaks are generally visible, and the smectic period d changes very little with P/L 
(1181 ). The reflectivity yields the electron density prohle along the director of the lamellar 


phase p{z) (averaged in-plane), but the analysis is rather involved (JIS 
shall not go into further detail here. 


27H30I) and we 


3.3 No interaction from a bilayer to the next 


A question of paramount importance is whether the pores interact from one bilayer to the 
next (along the z direction); we need to answer it in order to choose the theoretical model 
employed for describing the data (2D vs. 3D interaction) and, furthermore, to determine 
the biological relevance of our measurements. 


An effective way of determining the interbilayer interactions fjl3l ) is by measuring the 
scattering pattern at different swelling values and comparing the variation between 
different curves and against the expected form factor of the scattering object. We per¬ 
formed this investigation by exposing concentrated samples {P/L = 1/12.5) to different 
osmotic pressures, see section I^TTl Figure [5] shows detector scans through the peaks along 
(various symbols), as well as a cut through the simulated form factor in Figured] (red 
line). 

The hrst observation is that the measured curves are very similar; furthermore, their 
shape is qualitatively similar to that of the simulated form factor, if we neglect the presence 
of a slowly varying background, presumably due to thermal fluctuations of the lamellar 
phase (see next section). We can therefore conclude that there is no interaction between 
pores from one bilayer to the next. 

Thus, for the purpose of studying pore interaction, the bilayers in the solid-supported 
stacks we investigate can be considered independent, as one would require for modelling 
the cell membrane. Although the Ala monomers can be charged at neutral pH (M, Is^ . 
the 100 mM NaCl concentration (similar to that of biological media) reduces the Debye 
length to about 10 A, effectively screening the electrostatic potential; the only remaining 
interaction is that mediated by the bilayer. 
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3.4 Pore signal 

Figure [6] shows the detector scans along qy (out of the plane of incidence) for four different 
P/L values (indicated alongside the curves). A very intense and sharp component in 
gy = 0 (due to the specular beam) was removed for clarity. Scattering from the thermal 
fluctuations gives rise to a wide “bump” centered at the origin; to remove it, we fit the 
scans with a three-Lorentzian model (illustrated for the lower curve) and subtract the 
central component from the measured data. 

As a measure of fit quality we use the function divided by the number of points 
Apnts- The standard deviation for each experimental point is determined considering 
a Poisson distribution for the measured signal (before background subtraction) 


3.5 Hard disk model 


The simplest model for the interaction is that of hard disks confined in the plane. Using the 
“fundamental measure” approach, Rosenfeld (1331. Eq. (6.8)) provided a simple analytical 


expression for S'hd(gr )5 which is accurate over the entire concentration range we explore; 
the complete formula is given in the Appendix. 

First, the fits were performed for each scan individually, the hard disk radius R and 
the number density of pores n = 1/Apore being free parameters. For each scan, we tried 
the form factor for the 6-, 7- and 8-monomer pore. One fit example is displayed in Figure 
[7] (for P/L = 1/12.5), and the values of the fit parameters are shown in Figure [H]for all 


scans. 

The first conclusion is that the best agreement between the value of n obtained by 
fitting and that calculated using Equation ([T]) is obtained assuming a 7-monomer pore. 
The agreement is slightly worse for the hexamer and clearly off for the octamer; this can 
also be seen from the values of the function for the different individual fits (data not 
shown). We can therefore assume that we are dealing with 7-monomer pores. 

A very important result of the individual fits is that the value of R decreases with the 
P/L concentration from 24.8 to 17.9 A. One might understand an increase in radius at 
higher concentration due to the appearance of pores with more than seven monomers, but 
a decrease is clearly an unphysical result, which might indicate the presence of a “soft” 
repulsive interaction: as the concentration increases, the pores are forced closer together, 
overcoming this energy barrier. We therefore redid the fits including such a contribution. 

The samples with P/L = 1/12.5 at lower osmotic pressure yield sensibly higher values 
of R than that corresponding to c = 31% (see Figure El right). For reliability, we decided 
to ignore these points in further fits. This discrepancy does not correspond to a change 
in interaction between pores (see section E75]l : it originates most probably in the difficulty 
of obtaining a clear separation between the pore signal and the thermal scattering, which 
increases substantially with decreasing osmotic pressure (data not shown). 


3.6 Additional interaction 

We now consider a more complex interaction, consisting of hard core repulsion (when the 
pores are in contact) and an additional, longer-range term, corresponding to a bilayer- 
mediated interaction. For simplicity, we describe this component as a Gaussian and 
we account for its effect on the structure factor perturbatively, using the random phase 
approximation (RPA); see the Appendix for more details. 
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The experimental data are fitted simultaneously using the same parameters; R is the 
hard core radius, Uq corresponds to the amplitude of the additional component and ^ to 
its range (see Eq. (|5]) for the dehnition). 

Since the best individual hts were obtained with a 7-monomer pore, we impose the 
pore density calculated for this model (corresponding to m = 7 in Equation ([3])) as well 
as the form factor (blue curve in Figure [7]). We also checked that the £t quality (as given 
by the function) is better than for the 6- and 8- monomer pores. 

The £t results are shown in Figure [9l and the interaction potential is plotted in Figure 
[Tnl Comparison between the different £t configurations is detailed in Table [Hand Figure 
[721 Briefly, the presence of the additional interaction strongly increases the quality of the 
fit with respect to a £t with a fixed R but t/o = 0 (x^/-^pnts = 15.14, as opposed to 22.0)@. 
The fit quality is still much worse than that obtained by letting R vary with the P/L 
ratio, but in this latter situation more £t parameters are used, aside from the unphysical 
assumption of shrinking pore size. Even for varying R, the additional interaction yields a 
modest decrease in In this case we obtain a similar range ^ but a much lower amplitude 
t/o, most of the effect being “simulated” by the apparent R variation (see Table [H for the 
value of the htting parameters and Figure [TH for the plots). 

To summarize, we hnd that the interaction between 7-monomer pores of alamethicin in 
DMPC bilayers can be described by a hard core with radius 18.3 A, in excellent agreement 
with the geometrically estimated outer radius of the pore (18.2 A) and an additional 
repulsive interaction described by Eq. ([3]), with a range ^ = 31.5±0.27 A and an amplitude 
Uq = 4.74 ± 0.09 corresponding to a contact value U{2R) = 2AkBT. These very 
small error bars on the £t parameters should however be considered very carefully, since 
the most important source of error is probably the simplihed model for S{q). 


4 Discussion and Conclusion 


Very few experimental results point to the existence of lipid-mediated interaction between 
membrane inclusions; to our knowledge, they were all obtained by freeze-fracture electron 
microscopy (FFEM) and yielded directly the radial distribution function of the 


inclusions. The data were compared to liquid state theories fl38l-l40l) and could be described 


by a hard-core model with, in some cases, an additional repulsive or attractive interaction. 

In contrast, very sustained theoretical efforts aimed at understanding these systems 
started 30 years ago (41^44); they are either continuum-elasticity theories f 451-481 or 
more detailed models taking into account the molecular structure of the lipid bilayer 
()l4 Il5l. l4ll . l49l) . Two main origins for inclusion interaction have emerged, as discussed 
below. 


4.1 Hydrophobic mismatch 


A wide consensus has been reached as to the importance of “hydrophobic mismatch”, the 
difference in len gth between the hydrophobic part of the protein or peptide and that of the 
host membrane fl5nl.l5ll). However, the specihc way in which this mismatch is accomodated 


^Although this value seems very large, the fit quality is (visually) adequate and the error bars on the 
fit parameters quite small: the decimal places in Table [T] are significant. It is very likely that the Poisson 
distribution severely underestimates the standard deviation on each point. 





















for one particular system is not at all clear, especially when the peptide is longer than 
the lipid (the case of alamethicin in DMPC), since both bilayer compression/expansion 


and peptide tilt can be involved fl5lL [h 


The alamethicin/DMPC system was studied by NMR, hnding that the peptide is either 
parallel to the bilayer normal fl53l) or tilted by 10 — 20° fl54l) . conclusion supported by 
simulations (55). Moreover, the peptide exhibits a kink at the Pro^^ residue, (54, ^ 5^ . 
making the evaluation even more complicated. If one considers the entire pore as one 
(rigid) object, the tilt is probably very small, due to its size (j58|). Thus, the mismatch is 
likely compensated by bilayer expansion, which propagates over a few tens of A from the 
edge of the inclusion (58|, l59l) . values comparable to our experimental findings. 


4.2 Changes in lipid ordering 

Another -more subtle- effect is that an inclusion modihes the structure of the bilayer 
by perturbing the configuration of the lipid chains fll4 l4l[ l60l. roll) . In particular, the 


results of Lagiie et al. fll4 Il5l) are in semi-quantitative agreement with our observations: 
they extracted the lateral density-density response function of the hydrocarbon chains 
from the MD simulations of a DPPC bilayer fl62l) and used it to determine the interaction 
between “smooth” (no hydrophobic mismatch) hard cylinders embedded in the bilayer. 
For the largest cylinder radius they considered (9 A, about half that of alamethicin pores), 
they obtain a repulsive lipid-mediated interaction with a maximum value of 10 ksT and 
extending 20 A from contact (IJ). This study was followed by a comparison between 
different lipids, including DMPC (the lipid used in our experiments) f)l5l) : intriguingly, in 
this case they hnd a non-monotonic interaction, attractive close to contact and repulsive 
for larger distances. Furthermore, this interaction extends further than in the case of 
DPPC. We did not perform a more detailed comparison between their predictions and our 
experimental results, since the interaction potential varies considerably with the inclusion 
radius, but the agreement is certainly encouraging. 


4.3 Perspectives 

For a complete description of the perturbation and the interaction it induces, both hy- 
drophobic mismatch and changes in chain ordering must be taken into account (l63l. 16411. 


It has been pointed out repeatedly fl48l . ro4lj that the spontaneous curvature of the mono- 
layer radically changes the lipid-mediated interaction. To date, no consistent picture has 
emerged, due to theoretical difficulties but also to the lack of experimental data. 

The experimental work presented here consisted in determining the lipid-mediated 
interaction between alamethicin pores in DMPC bilayers; we found it to be repulsive 
and the overall shape of the potential is in qualitative agreement with recent theoretical 
predictions (0,0. However, the quality of the fits to the experimental data is not very 
good; this can stem from technical difficulties and systematic errors, but also from the 
rough model employed (RPA approximation). Both these aspects will be improved in the 
future but the results are already significant. 
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Appendix 


Hard disk model 


We nsed the analytical expression for the strnctnre factor of hard disks given by Rosenfeld 
a, Eq. (6.8)): 


>5hd(?) 


1 + 4?7 


(Ji{qR)y , ^MqR)Ji{qR) , ^Ji(2gR) 


(4) 


where q is the in-plane scattering vector, R the hard disk radins, p = mrR^ the packing 
fraction (with n the nnmber density of the disks) and Jk the Bessel fnnctions fnnctions 
of the first kind and order k. The prefactors are given by: 


G = (l-r/)-3/2 
_ l + r] 

^ - (TV^ 

A = r]~^ [1 {2r] - l)x + 2r]G] 

B = r]-^ [(1 -r])x-l- 3r]G] 

Additional repulsive interaction 

We added a repnlsive component described by a Ganssian, with amplitnde Uq and range 

e: 


U{r) 


Uq exp 



(5) 


considered as a pertnrbation with respect to the hard disk model, taken into acconnt via 
the random phase approximation (RPA) fjosl) . In this approach, one obtains the direct 
correlation fnnction of the pertnrbed system c(r) from that of the reference system Cref(r) 
as 


(1661 ) or, eqnivalently: 


c(r) = Cref(r) - /3U{r) 


S-\q) = S;J{q)+pf3U{q) 


with U{q) 


2ttUq exp 


M! 

2 


the Fonrier transform of f/(r). 


( 6 ) 

(7) 
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Fit parameters 

For all practical purposes, we give in Table [1] the value of the £t parameters for to the 
different conhgurations discussed in the text; throughout, the form factor and the density 
are those of a 7-monomer pore. The corresponding hts are displayed in Figure [121 The 
top left set (same R, same Uq) is the same as in Figure 0 
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Tables 


Table 1: Fit results with different models. Fit conditions refer to the hard core radius R being 
the same for all scans or allowed to take different values for different individual scans (‘free’) and 
to the presence or absence (‘Co = O’) of the additional interaction. Uq and ^ are the amplitude 
and range of the additional interaction. x^/Apnts is an indication of the fit quality. Aparam is 
the number of fit parameters, including the seven intensity prefactors, one for each scan. 


P/L 

Par am. 

Same R 
same Uq 

Same R 
Uo = 0 

Free R 
same Uq 

Free R 
Uo = 0 

1/25 

R[A] 

18.3 

19.2 

24.3 

24.8 

1/20 

R 

D 

D 

23.5 

24.1 

1/15 
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D 
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21.7 

22.1 

1/15 
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D 
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22.8 

23.1 

1/12.5 
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D 

D 

21.4 

21.7 

1/10 
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D 
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19.4 

19.8 

1/7.5 
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D 

D 

17.7 

17.9 


Uo [kBT] 

4.74 

0 

1.56 

0 


e[A] 

31.5 

- 

34.8 
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N 

’ param 

10 

8 

16 

14 


xV^pnts 

15.14 

22.0 

8.93 

9.44 
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Figures 
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Figure 1: Absolute value of the form factor \F{qy,qz)\ for a 7-monomer alamethicin pore after 
subtraction of the pure bilayer background and azimuthal averaging. 
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Figure 2: Reconstruction of the electron density profile of the 7-monomer alamethicin pore 
(after subtracting the pure bilayer background) by Fourier transforming back to real space the 
form factor displayed in Figured) 
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Figure 3: Structure of the reciprocal space for a lamellar stack. Reflectivity scans are performed 
along the vertical z axis, while the CCD images -an example of which is shown to the right- 
represent slices through the reciprocal space along the Ewald sphere (shown in red). The charac¬ 
teristic featnres can easily be identified: the intense diffuse sheets around the very sharp Bragg 
peaks; the extended and narrow circle arcs going through the Bragg peaks are defect-induced 
Debye-Scherrer rings. The intensity increase close to the horizon is due to dynamic effects. 
Finally, the off-axis signal exhibiting a maximum at about qy = 0.1 A is due to the presence 
of alamethicin pores. 



Figure 4: Reflectivity spectra of aligned DMPC multilayers containing alamethicin. The data 
is only shown for four P/L concentrations. Curves vertically shifted for clarity, with the P/L 
ratio increasing from bottom to top. 
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Figure 5: Sections along qz through the pore signal of samples with P/L = 1/12.5 at different 
osmotic pressures (symbols) for qy = 0.1 A ^ and cut through the square of the simulated form 
factor \F{qz)\‘^ (red line). The sharp peaks appearing in the top and bottom curves are due to 
the Debye-Scherrer rings (see Figure [3]) 



Figure 6: Detector scans along qy through the pore signal for four P/L concentrations (in¬ 
tegrated in the range 0.14 < ( 7 ;^ < 0.18 A ). For the bottom scan we also show the three- 
Lorentzian fit to the data; the central “bump” is subtracted before further treatment. 
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Figure 7: Fit of the data for P/L = 1/12.5 with a hard-disk model. Diamonds: point detector 
scan. The fit, shown as a red line, is the product of the form factor \F{qz)^ for the 7-monomer 
pore (blue line) the structure factor for a hard disk system (black dashed line), with radius 
R = 2 O. 95 A and number density n = 1/^pore = 3.59 10“^A . For clarity, only the qy > 0 range 
is displayed. 
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Figure 8: Values of the fit parameters (density and hard disk radius) obtained from individual 
fits of the scans, with the form factor of 6- (open diamonds), 7- (solid dots) and 8-monomer pores 
(open triangles). For comparison, the density obtained as n = 1/Apore according to formula ([I|) 
is shown as dotted, solid and dashed line, for the 6-, 7- and 8-monomer pore, respectively. 
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Figure 9: Experimental data (dots) and fits (solid lines) with a hard disk model and an addi¬ 
tional repulsive contribution. The curves correspond to different peptide concentrations; from 
top to bottom, Ala/DMPC= 1/7.5, 1/15, 1/25, 1/20, 1/12.5, 1/10, 1/15. All scans are fitted 
simultaneously, yielding R = 18.3A, Uq = A.lAksT and ^ = 31.5 A. 



Figure 10: The interaction potential used for the fits in Figure [9l consisting of a hard core and 
an additional Gaussian repulsion, given by formula ([S]), with parameters Uq = and 

^ = 31.5 A. The contact value U{2R) = 2AlkBT. 
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Figure 11: Schematic representation of interacting pores in a lipid bilayer. The seven monomers 
(red) border the water pore (blue). In grey, the range of the lipid-mediated repulsion. 
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Fit 1: same R, same Uq 



Fit 2: same R, Uq = 0 
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Fit 3: free R, same U, 
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Figure 12: Experimental data (dots) and fits (solid lines) in the various configurations described 
in Tabled! The top left set, corresponding to the same R and Uq, is the same as in Figure [9| 
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